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Abstract 

We systematically study the most general Lorentz-violating graviton mass invari- 
ant under three-dimensional Eucledian group. We find that at general values of mass 
parameters the massive graviton has six propagating degrees of freedom, and some of 
them are ghosts or lead to rapid classical instabilities. However, there is a number 
of different regions in the mass parameter space where massive gravity is described 
by a consistent low-energy effective theory with cutoff ~ \JmMpi. This theory is free 
of rapid instabilities and vDVZ discontinuity. Each of these regions is characterized 
by certain fine-tuning relations between mass parameters, generalizing the Fierz-Pauli 
condition. In some cases the required fine-tunings are consequences of the existence 
of the subgroups of the diffeomorphism group that are left unbroken by the graviton 
mass. We found two new cases, when the resulting theories have a property of UV 
insensitivity, i.e. remain well behaved after inclusion of arbitrary higher dimension 
operators without assuming any fine-tunings among the coefficients of these operators, 
besides those enforced by the symmetries. These theories can be thought of as gener- 
alizations of the ghost condensate model with a smaller residual symmetry group. We 
briefly discuss what kind of cosmology can one expect in massive gravity and argue 
that the allowed values of the graviton mass may be quite large, affecting growth of pri- 
mordial perturbations, structure formation and, perhaps, enhancing the backreaction 
of inhomogeneities on the expansion rate of the Universe. 



1 Introduction 

During the last few years a significant progress in understanding of properties of massive 
gravity in four dimensions has been achieved. One of the first steps was the implementation 
of the Stiickelberg trick in the Fierz-Pauli theory [1]. As a result, a clear unified picture 
of the apparently mysterious properties of the massive gravity has emerged. Namely, the 
explicitly covariant Stiickelberg language made it clear that the famous van Dam-Veltman- 
Zakharov (vDVZ) discontinuity [2, 3] and related non-linear effects [4], as well as the hidden 
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sixth polarization mode of the massive graviton [5] are all consequences of the fact that the 
kinetic term in the Goldstone sector of the Fierz-Pauli theory is degenerate in the limit 
Mpi — > oo. As a result, the theory is not well defined in the limit when gravity decouples 
and the cutoff energy of the Fierz-Pauli theory in the Minkowski background is extremely 
low, 

As < (M 2 Pl m) 1/z ~ (1000 km)" 1 , 

where m is a graviton mass, which we take of order the current Hubble scale. It is worth 
noting that this may not necessarily imply that the Fierz-Pauli theory is unacceptable phe- 
nomenologicaly. Indeed, recently it was shown [6] (see also [7]) how similar problem [8, 9] 
present in the context of co-dimension one Dvali-Gabadadze-Porrati (DGP) model [10] may 
be cured by taking into account the effects of local curvature. It is unclear at the moment 
whether this mechanism can be implemented in the Fierz-Pauli case. 

The further progress in the four-dimensional theories of massive gravity is due to idea [11], 
that instead of Fierz-Pauli mass term, which completely breaks diffeomorphism invariance 
of the theory, one can introduce a graviton mass, breaking only the time reparametrization 
invariance, 

t^t + ?(t,x) , 

while keeping the arbitrary time- dependent spatial diffeomorphisms 

x'^x' + C^x'). (1) 

The crucial difference with the Fierz-Pauli case is that in this model the Goldstone sector 
(consisting of a single scalar field) is well behaved in the limit when gravity is decoupled. At 
the linear level the resulting Lorentz non- invariant theory (dubbed "ghost condensate") is 
characterized by two parameters — the characteristic length and time scales. Interestingly, it 
appeares that both of these parameters can be substantially shorter than the inverse Hubble 
scale without contradicting to observations [12, 13] 1 . It is worth mentioning, however, that 
in many respects ghost condensate is very different from what one would call a theory of 
massive graviton. In particular, in this model a massless transverse traceless mode is present 
in the spectrum. 

1 It was argued in Ref. [14] that the limits can be much tighter due to efficient accretion of the ghost 
condensate by black holes. It is not clear at the moment, whether the steady state solutions found in 
Ref. [14] are the physically relevant ones. 
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Finally, very recently it was argued in Ref. [15], that the generic Lorentz- violating gravi- 
ton mass term, preserving Euclidean symmetry of the three-dimensional space 2 , 

M 2 

L m = {mlh 00 hoo + 2m\h 0i h 0i - m 2 /^/^- + mlhahjj - 2m\h 00 hn) , (2) 

where are perturbations about Minkowski metric, also leads to the theory with well 
behaved Goldstone sector in the decoupling limit. More precisely, it was shown to be true 
for m = and provided the other masses satisfy certain positivity conditions. As a result, 
this mass term does not lead neither to ghosts, nor to low strong coupling scale in the flat 
space. Unlike the ghost condensate model, the spectrum of the theory has a mass gap, and 
describes five propagating degrees of freedom. Still, the vDVZ discontinuity is absent in this 
model and the cutoff scale is given by 

A 2 ~ y/mMpi ~ (0.1mm)" 1 . (3) 

However, one may be worried by the need to satify the fine-tuning relation 

m 2 = 

in this theory. This condition is remininscent of the Fierz-Pauli fine-tuning and may be a 
signal of the presence of the analogue of the Boulware-Deser mode and related instabilities. 

This result indicates that the ghost condensate is unlikely to be a single consistent infrared 
modification of gravity and calls for the further systematic study of the Lorentz-violating 
graviton mass. This is the main task of the current paper. 

We start in section 2 with the covariant formulation of the Lorentz-violating massive 
gravity. For this purpose we use a Stiickelberg-like formalism. This formalism is a slightly 
more flexible version of one used in the Lorentz-invariant case [1]. As in the Lorentz-invariant 
case, graviton mass generically leads to the emergence of four new dynamical degrees of free- 
dom — Golsdtone bosons corresponding to broken reparametrizations of four coordinates. 
An example of ghost condensate hints that the Goldstone sector can be made well behaved 
in the limit when gravity is decoupled by requiring that some residual subgroup of the diffeo- 
morphism group is left. However, the group of diffeomorphisms has too many subgroups, so 
the case by case consideration of all possible patterns of symmetry breaking seems hopeless. 

Instead, in section 3 we adopt the bottom-up approach and study the spectrum of the 
Lorentz-violating massive gravity at general values of masses in the decoupling limit 
2 A possibility to introduce mass terms of this type was also mentioned in Ref. [16]. 
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Mpi — > oo. In agreement with expectations of Ref. [15] we find that at non-zero itiq theory 
generically has either ghosts or exhibits rapid classical instabilities. However, we found a 
number of other fine-tuning relations ensuring the absence of pathologies in the Goldstone 
sector at the one-derivative level. Actually, the very possibility of the other choices was 
expected, as we know that region with > and all other masses equal to zero correspond 
to the ghost condensate model with a well behaved Goldstone sector. What is interesting is 
that one obtaines the well behaved Goldstone sector under weaker assumptions as well. 

In section 4 we discuss the effects of higher derivative terms. We argue, that, generically, 
one should assume an infinite number of fine-tunings in higher order terms to preserve the 
above attractive properties of the low-energy theory. It is worth stressing that this does not 
imply that these fine-tunings cannot be realized in some specific UV completions of massive 
gravity. This only means that a generic UV completion would result in a massive gravity 
theory with cutoff scale much lower than A2. 

However, we found two new cases when all the required fine-tunings can be ensured by 
the residual reparametrization symmetry. In both cases this symmetry is a subgroup of the 
residual symmetry present in the ghost condensate. Namely, in one case massive gravity is 
invariant under time independent spatial diffeomorphisms 

and in other under time-dependent shifts of spatial coordinates 

X ^^ + C(t). (4) 

Both theories inherit from the ghost condensate the property of UV insensitivity. Namely, 
they remain well behaved under inclusion of all possible higher derivative operators sup- 
pressed by the cutoff A 2 and compatible with the residual reparametrization symmetry, 
without assuming any fine-tuning relations between their coefficients. 

Some of other properties of these massive gravities are also similar to that of the ghost 
condensate. For instance, all three models have a scalar degree of freedom with peculiar 
dispersion relation 

2 4 

CO (X p . 

However, in many respects two new models are very different from the ghost condensate. 
For instance, the tensor graviton mode has non-zero mass in case of residual symmetry (4). 

In section 5 we discuss effects due to mixing with gravity, and prove that vDVZ dis- 
continuity is absent in the massive gravity models discussed here (both fine-tuned and UV 
insensitive). 
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We conclude in section 6 with a brief discussion of cosmology in massive gravity. We pro- 
vide an example of a theory with massive tensor mode, where flat homogeneous cosmological 
solutions are identical to those in the Einstein theory. This implies that the graviton mass 
can be substantially larger than the current Hubble value. Consequently, graviton mass may 
modify in the interesting way the dynamics of primordial perturbations and structure forma- 
tion or, even, gravitational dynamics at the (super)galactic scales. Even more intriguingly, 
there is a possibility, that the presence of small scale inhomogeneities may significantly affect 
the expansion rate in massive gravity, potentially providing a link between the onset of the 
cosmological acceleration and epoch of structure formation. 



2 Stiickelberg trick for Lorentz- violating massive grav- 
ity 



In the most general sense by massive gravity we understand any theory described by the 
following action 



where the first term is a usual Einstein-Hilbert term and F is, generally speaking, an ar- 
bitrary function of metric components, their derivatives, and coordinates t, x l (i = 1, 2, 3) 
itself. One can add some matter fields to this system, which we assume to be minimally 
coupled to the metric. In the current paper we will limit our consideration to the case 
when theory (5) admits empty Minkowski space-time as one of the solutions to the field 
equations. Also, following Ref. [15], we assume that a low-energy effective theory for small 
perturbations around this solution is invariant under the Eucledian transformations of the 
three-dimensional space coordinates x l . In particular, this implies that function F has no 
explicit dependence on the coordinates. Furthemore, we assume that function F depends on 
a single energy scale A. 

Then, the quadratic theory near Minkowski background is described by the usual lin- 
earized Einstein-Hilbert action plus mass term of the form (2), where the overall mass scale 
m is related to the parameters A and Mp\ as 



The key issue we address in the current paper is what conditions function F should satisfy 
in order dynamics of metric perturbations near Minkowski background could be described 




(5) 



A = ^mMpi . 
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by a low-energy effective field theory with cutoff A. It was shown in Ref. [1], that this is 
not possible in the Lorentz-invariant case. The convenient tool allowing to adress this issue 
is the Stiickelberg formalism. Our version of this formalism is essentially equivalent to that 
used in Ref. [1]. Let us nevertheless briefly describe it, both to make our presentation more 
self-contained and because we are using slightly more flexible version of the formalism to 
allow for Lorentz non-invariant situation. 

The first step is to introduce four scalar fields M (/i = 0, ... ,3). These are Goldstone 
bosons corresponding to four broken diffeomorphisms 

The dynamics of these fields is described by some sigma-model action. We assume that this 
action is characterized by a single energy scale A. The main requirement to this action is 
that when coupled to gravity it provides the following field configuration as a solution to the 
field equations 

^ = A V (6) 
9tMu = G^ v , (7) 

where G^ u is some fixed metric (Minkowski metric in our case). Provided such a solution 
exists, one can write 

9txv = + (8) 

= A V + 7T^ (9) 

and study dynamics of small scalar and metric perturbations 7r M , h^. This, by definition, 
will be a certain version of massive gravity in the space-time with metric G^ v . Indeed, given 
such a theory, one can fix the unitary gauge 7r M = and recover the previous non-covariant 
formulation (5) of massive gravity 3 . 

Inversely, given action (5) one always can make it explicitely covariant reintroducing the 
Goldstone fields M . Namely, everywhere in the action (5) one replaces metric g^ v with its 
gauge transform 

g^x) - 9» v ee Q y d y 9Xp (Y(x)) . (10) 
Here are component functions of the diffeomorphism 

Y : x 1 * -> Y^x) . 

3 Note, that field equations following from the variations of the Goldstone fields are consequences of the 
other equations in the unitary gauge. 
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Analogously, in the presence of matter fields ip one replaces them everywhere with their 
images under gauge transformation. Clearly fields appear only in the non-covariant 
part of the action. In particular, no direct coupling between these fields and matter arises, 
provided the action for matter is covariant. 

Then action (5) considered as a functional of the metric g^ u and fields is diffeo- 
morphism invariant under conventional transformations of the metric, accompanied by the 
following transformations of the fields Y^, 

Y^ir'otr, (ii) 

where £ is a change of coordinates and o is a natural multiplication of two diffeomorphisms. 
Transformation law (11) implies that component functions do not transform as scalars 
under changes of coordinates. However, the transformation law for the component functions 
of the inverse mapping Y~ l is that of the scalar fields, 

(y-r ■ 

These component functions are precisely the Goldstone fields 0^ introduced earlier, 

Now, imagine that massive gravity in the unitary gauge (5) has a residual invariance un- 
der diffeomorphisms from a certain subgroup 7i of the whole group of diffeomorphisms. 
In covariant formalism this implies that the theory is invariant under the following transfor- 
mations of component fields Y^, 

Y^{Yoi H y. (12) 

In terms of scalar Goldstone fields M this implies that a group Q of internal global symmetries 
of the Goldstone action contains group Ti as its subgroup, i.e. transformations of the type 

r - e H {<t>) (is) 

are symmetries of the Goldstone action. One could think that it is impossible to con- 
struct non-trivial actions with such symmetries for large enough groups TC. However, the 
Stiickelberg trick described above allows to construct a lot of such actions starting from the 
massive gravity in the unitary gauge. 

If, in addition, a group of internal symmetries Q contains the common subgroup with the 
group of isometries of the metric G^ v , then vacuum (6), (7) is invariant under the diagonal 
combination of these two subgroups. 



7 



To finish this general discussion the following comment is in order. One may question the 
assumption that the massive gravity contains just metric and Goldstone bosons. Indeed, in 
the case of gauge fields Higgs mechanism implies the presence of at least one additional scalar 
field — the radial mode of the Higgs field. Presumably, the successful implementation of the 
Higgs mechanism in gravity (i.e., a theory of massive graviton with UV cutoff at the Planck 
scale) would also require some extra fields. Unfortunately, such a theory is not constructed 
yet 4 . In this situation, the best we can do is to introduce the minimal set of fields which is 
required by the consistency. Then there is nothing wrong in solving the resulting non-linear 
equations, provided one restricts oneself to the classical solutions tractable within the low- 
energy effective theory. Tractable means here, that the theory of small perturbations around 
the solution has no pathologies leading to rapid instabilities and has large enough cutoff 
scale. Solutions that do not satisfy these criteria cannot be studied without more detailed 
knowledge of the UV completed theory. 

Let us illustrate the above general discussion by some concrete examples. For instance, 
the Fierz-Pauli theory corresponds to the choice of the Minkowski metric as a background 
metric G^ u , and the four-dimensional Poincare group (i.e., the whole group of isometries of 
the Minkowski space-time) as a group 5 Q. 

Ghost condensate corresponds to the degenerate situation when action does not depend 
on the fields 0* at all. Formally, one can say that in this case a group of internal symmetries 
Q is a subgroup of the diffeomorphism group shifting time by a constant and transforming 
spatial coordinates in the arbitrary way. 

The Lorentz- violating massive gravity considered in Ref. [15] is also a theory of massive 
gravity in the Minkowski space-time. It corresponds to the choice of the subgroup of the 4d 
Poincare group generated by space-time shifts and rotations of space coordinates as a group 
of internal symmetries. Then the most general Goldstone action has the following form at 
the one-derivative level 



4 See, e.g., Ref. [17] for a recent attempt in this direction, making use of the vector field. 

5 Strictly speaking, in this way one obtains some general Lorentz-invariant theory of massive graviton. To 
ensure that the graviton mass has the Fierz-Pauli form, the sigma-model action should satisfy additional 
fine-tuning conditions. 




(14) 
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where X, F lJI , V 1 and Z are the following scalar quantities 

X = fTd^d^/A 4 
Y ij = g^d^d^/A 4 
V i = g^d^d^/A 4 
Z = e^d^d^dxtfd^/i^A 8 ) (15) 

In the rest of the paper by massive gravity (or Lorentz-violating massive gravity) we imply 
this particular model. Note that Z 2 can be written as a function of X, Y l \ V 1 , so that 
different functions F in Eq. (14) may lead to the same Goldstone action. The Latin indices 
i,j = 1, 2, 3 are converted using the 3d Kronecker symbol 5ij in Eq. (14). The energy scale 
A is a naive UV cutoff of this model which is assumed to be much lower than Mp\. As we 
explain later, in certain circumstances the actual UV cutoff may be significantly lower. 

Clearly, for constant metrics g^ v the field configuration (6) is always a solution to the field 
equations following from the action (14). To satisfy the second condition (7) for a Minkowski 
space-time one should find two positive constants a, b > such that the energy-momentum 
tensor of the sigma-model (14) is zero in the background metric 

r)^ = diag(a, -b, -b, -b) . 

This condition translates into a following pair of equations 

2F X X = \F Yii Y ij = F-ZF Z (16) 

for the values of function F and its first derivatives 

dF _ dF dF 

Fx = dx' Fy " = W^> Fz = dz 

at the point 

X = a' 1 , Y l] = -b- 1 ^ , V 1 = , Z = a~ 1/2 b~ 3/2 . 

It is straightforward to check that for a generic function F one can always find a and b such 
that these conditions are met. Then masses m^, . . . ,m\ are combinations of the first and 
second derivatives of the function F at this point. The overall mass scale m is related to A 

as 

A ~ yjMpim . 
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3 Quadratic Goldstone action at the one-derivative level 



Let us first study the Goldsone sector for general values of masses m\ in the absence of 
gravity. To obtain the quadratic Goldstone Lagrangian one plugs the "pure gauge" metric 
perturbation 

in the action (2). The result is [15] 

L = M 2 Pl (2m 2 (d iT ) 2 + m 2 (doiri) 2 + m 2 (<9j7r ) 2 + (4m 2 - 2ml)ir d d i ir i - m 2 ^^-) 2 - 

(m 2 2 -2mj)(d z 7c z ) 2 ) ■ (17) 

This Lagrangian can be partially diagonalised by decomposing fields into components cor- 
responding to different irreducible representations of the three-dimensional Eucledian sym- 
metry group. The vector sector contains transverse part nj of the fields n\ 

d\ T = . 

The high-energy action in this sector depends only on masses m\ and m| and has the 
following form 

L = M 2 Pl (m 2 1 (d nf) 2 -m 2 2 (d l ^) 2 ) ■ (18) 
In general, this action is well behaved for 

mj, m\ > . (19) 

Indeed, if ml is negative, then nf become ghosts, while if ml is positive and m\ is negative 
then there are classical instabilities at all spatial momenta. Points where at least one of 
these masses is zero are special. Thus conditions 

ml = (20) 
m\ = (21) 

are candidates on the role of fine-tuning relations. At the moment we proceed under assump- 
tion that condition (19) is satisfied and discuss later what happens if one or both of Eqs. (20), 
(21) hold. It is worth recalling [15] also, that m\ is a mass of the transverse-traceless tensor 
graviton mode. 

The scalar perturbations are 7r and ttl, where the latter is defined as 

7Tj = djlYL ■ 

V^d 2 
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Their Lagrangian has the following form 



L = M 2 Pl {2ml{d 7r ) 2 + ml{d 7r L ) 2 + m^7r ) 2 - (Ami ~ 2ml)n d ^ -d 2 n L - 

2(mi-mj)(d t n L ) 2 ) . (22) 

In the Fourier space the field equations following from this Lagrangian are 

Ml n ° 1=0, (23) 
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where M is the following 2x2 matrix, 



/ au 2 + 1 —iiiv i 
M =[ ■ 2 R . (24) 



where 

a; 

z/ = — 
P 



2m^ 2ml — ml R m 2 ' '" 

Q; = 2~ j A 4 = 2 ' I = 2 
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m\ m\ 



with c<j and p being energy and absolute value of the three-momentum correspondingly. 
Generically, there are two propagating degrees of freedom in the scalar sector. Their disper- 
sion relations are determined by the condition that the determinant d of the matrix M is 
zero, 

d = au A + u 2 (l-Pa-fx 2 )-p = 0. (25) 

There are some special values of the parameters, where the number of degrees of freedom in 
the scalar sector is reduced. It may happen either when one of the roots of Eq. (25) goes to 
infinity, so that the dispersion relation for the corresponding mode takes form 

P 2 = , (26) 

or when one of the roots goes to zero, so that the corresponding dispersion relation is 

u 2 = . (27) 

The first option is realized for 

a = (28) 
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so that the order of the equation (25) reduces by one. This is the case considered in Ref. [15]. 
There is a single propagating mode in the scalar sector in this case. Actually, there is a 
possibility to reduce the order of the on-shell equation even further. This happens when 

a = , /i 2 = 1 . (29) 

In this case there are no propagating degrees of freedom in the scalar sector. The second 
option (27) is realized for 

P = 0. (30) 

Both scalar degrees of freedom have zero energy and do not propagate if 

= , /i 2 = 1 . (31) 

Finally, it may happen that one of the fine-tuning relations (28), (29) and one of (30), 
(31) are satisfied simultaniously. For instance, the Fierz-Pauli fine-tuning implies that both 
conditions (29) and (31) hold (with /j, — 1). 

Actually, we are making things somewhat too simple saying that the above special cases 
correspond to the reduction in the number of degrees of freedom. This is literally the case 
at the quadratic level in the limit when gravity is decoupled and in the absence of higher 
derivative terms in the Goldstone sector. However, one may be worried that inclusion of any 
of these effects strongly modifies the whole picture. This really happens in some cases, and 
more detailed discussion of these effects will be presented in sections 4, 5. The oversimplified 
picture is enough, however, for the main purpose of this section — to identify regions in the 
mass space, where obvious pathologies are absent in the low-energy effective theory. Let us 
now perform case by case analysis of different possibilities. 



3.1 Generic graviton masses: six propagating degrees of freedom 
and ghosts or classical instabilities 

To start with, let us consider the most general case when none of the above fine-tuning 
relations holds, so that there are two propagating degrees of freedom in the scalar sector, and 
massive graviton has six polarization modes in total. Let us check whether it is possible then 
to chose parameters a, /x, f3 in such a way that there are no ghosts or classical instabilities. 

Formally the problem is formulated as follows. At any given values of the parameters a, 
/i, P and positive v 2 the Hermitian matrix M has two real eigenvalues A±(z/ 2 ) given by 
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where t and d are trace and determinant of the matrix M. Also there are two values v\ 2 
where the on-shell condition (25) holds, i.e. one of the eigenvalues is zero at this point. 

The question is whether it is possible to find values of the parameters such that the 
following two conditions are met, 

• Absence of the classical instability, 

v\ 2 > (33) 

• Absence of ghosts 



dv 2 



> (34) 



Let us prove that this is impossible. The first observation is that to avoid ghosts it is 
necessary that each of the eigenvalues A±(i/) has one zero. Indeed, if one of them has two 
zeros and another none than the derivative in one of the zeros is necessarily negative. Similar 
situation happens in a theory of a single scalar field with four derivative kinetic term. 

Now, note that Eq. (25) implies that one of the conditions for the absence of the classical 
instability is 

^ < (35) 
a 

Let us first discuss the case a < 0, (3 > 0. Here, in the limit of large v 2 the determinant of 
the matrix M is 

d = au 4 + 0{y 2 ) < 

Consequently, the smaller eigenvalue A_ is negative at large v 2 in this case. This implies, 
that derivative of A_ in its zero is negative, and this zero corresponds to the ghost. 
In the another case a > 0, (3 < the trace of the matrix M, 

t = (a + \)v 2 + 1 - (3 (36) 

is positive definite. This implies, that A + is also positive definite. Consequently, A_ has two 
zeroes, and one of them corresponds to the ghost. 

3.2 Phase m 2 = 

Let us now consider the case when the only fine-tuning relation which holds is (28). This is 
the case considered in Ref. [15]. 
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For a = one of the solutions v 2 to the on-shell condition (25) becomes infinite. The 
only remaining solution is 

JL 

l-fi 2 



-o = ~t~2 (37) 



Let us first consider the case 



Then posit iveness of z/q implies that 



/i 2 > 1 . 

(3 < . 



As before, in this case the trace (36) of the matrix M is positive for all positive z/ 2 , so A + is 
also positive and z/q is a zero of A_. On the other hand the determinant d of matrix M is 
now equal to 

d = i/ 2 (l - /.I 2 ) - (3 (38) 

and negative at large v 2 . Consequently, A_ is negative at large v 2 and has a negative 
derivative in its zero z/q. So the only dynamical degree of freedom in the scalar sector is a 
ghost here. 

Let us now consider the opposite case 

p? < 1 . (39) 

Then in order to avoid classical instabilities (3 should be positive, 

(3 > . (40) 

In this case both determinant (38) and trace (36) are positive at large z/ 2 . Consequently, 
both eigenvalues A ± are positive there as well. This implies, that z/ 2 is a zero of the smaller 
eigenvalue A_, and that derivative of this eigenvalue in its zero is positive. 

Furthemore, if in addition to Eq. (28) one also has /i 2 — 1, then the second root of 
the dispersion relation (25) goes to infinity. The determinant of matrix M is equal to —(3 
and non-zero in this case, so one can use the inverse of this matrix as a propagator in the 
perturbation theory (for (3^0). The propagating scalar degrees of freedom and obvious 
pathologies are absent in the Goldstone sector in this case. If (3 is equal to zero, then both 
fine-tuning relations (29) and (31) hold and we are in a situation generalizing that in the 
Fierz-Pauli theory. Here the quadratic part of the Goldstone sector is zero, so one cannot 
organize perturbation theory in this sector in the limit when gravity is decoupled. 

Finally, there is a possibility that fine-tuning conditions (28) and (30) are satisfied simul- 
taneously. Here the second mode in the scalar sector has dispersion relation of the form (27) 
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and is not propagating as well. There are no obvious pathologies in the Goldstone sector in 
this case, provided /j 2 7^ 1. 

To summarize, Goldstone sector is free of obvious pathologies in the high-energy limit at 
m l — 0) provided the graviton masses satisfy the following conditions 

ml = 0, mf, ml > 0, m 2 > m 2 > 0, m 2 > m 2 

or 

ml = 0, ml, ml > 0, (m 2 — m%)m\ = 0, ml 7^ ml 

or 

ml = 0, ml, ml > 0, m 2 = m 2 , (m 2 — m 2 )m 2 7^ . (41) 

In the bulk of this phase one of the modes in the scalar sector has dispersion relation (26) and 
is not dynamical, so in total there are three propagating degrees of freedom in the Goldstone 
sector. 

We will discuss later what happens when either m 2 or ml is zero. Already now it is clear 
that if in any of the above situations one has 

ml = , 

then nothing changes in the scalar sector, while both vector degrees of freedom become 
non-dynamical and tensor mode is massless. 

The above conditions are in agreement with those found in Ref. [15]. The difference is 
that we also studied here what happens on the boundaries of the region found in that work. 
Another difference is that an extra condition Ami > m l was imposed there. We will discuss 
the origin of this extra condition and whether it is really necessary in section 5. 

3.3 Phase ml = ml 

In this subsection we proceed under assumption that both conditions (19) are satisfied, and 
study what happens if ml is non-zero and fine-tuning relation (30) holds. Determinant d of 
the matrix M takes the following form in this case 

d = u 2 (au 2 + 1- /1 2 ) (42) 

while its trace t is 

t = (a+ \)v 2 + 1 . (43) 
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One of the solutions of the on-shell equation (25) is v 2 = and does not propagate. Note, 
that this mode always corresponds to a zero of the smaller eigenvalue A_, because trace (43) 
is positive at v 2 = 0. The second solution of the on-shell equation is 



Now, for a < the determinant (42) is negative at large z/ 2 , and, consequently, A_ is also 
negative there, while A + is positive. The only possibility to have v\ z positive in this case is if 
z/| 3 is also a zero of A_, and if the derivative of A_ is negative at u 2 3 . Hence, one inevitably 
has classical instability or ghost for m\ = ml and m\ < (assuming conditions (19) are 
satisfied). 

For a > both trace (43) and determinant (42) are positive at large v 2 , and consequently 
A + and A_ are also both positive there. Therefore, provided z/| 3 is positive, it is a zero of 
A_ in this case, the sign of the corresponding derivative is positive and there are no obvious 
pathologies in the Goldstone sector. 

We have already discussed what happens if a — 0. Another possibility is that /i 2 = 1, 
i.e. fine-tuning relation (31) is realized. In this case, there are also no obvious pathologies 
in the Goldstone sector. To summarize, there is a phase 

m\ = m 2 > , mf > , m\ > , m 2 > m\ or m\ < 

or 

ml = ml > , m 2 > , (m\ - m 2 )m^ = , m\ ^ , (45) 

where Goldstone sector is well defined in the limit when gravity decouples. In the bulk of this 
phase one of the modes in the scalar sector has dispersion relation (27) and is not dynamical, 
so in total there are three propagating degrees of freedom in the Goldstone sector. 

As above, nothing changes in the scalar sector if one takes in additon m\ = so that 
there are no propagating degrees of freedom in the vector sector and tensor mode is massless. 

3.4 Phase ml = 

Let us discuss now what happens when mass m\ = 0, so that gradient terms for the vector 
modes vanishes and tensor modes are massless. The above discussion has already covered 
this case for positive m\. If, in addition to m|, m\ is also zero, then the whole quadratic 
Lagrangian in the vector sector is zero in the limit when gravity decouples, so the Goldstone 
sector is not well behaved in this limit (the only exception is the ghost condensate case, 
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when there are unbroken time-dependent spatial diffeomorphisms, so that Goldstones it 1 do 
not show up at all). So the only remaining possibility is that m\ is negative, 

m\ < (46) 

This is possible because vector modes have zero energy and do not propagate for m\ = 0. 

The analysis of the scalar sector proceeds now very similar to what has been done above. 
The only difference is that the condition that there are no ghosts in the scalar sector is now 
opposite to (34), 



dv 2 



< . (47) 



We will skip technical details and present just the results. As above, there are pathologies 
in the scalar sector for generic values of the parameters a, //, f3. For a = the allowed region 
is 

or 

m o = m 2 = ? m i < ? m 3 = , 777.4(777.4 — 777^) 7^ 

or 

m o = m 2 = > ml < , m\(m\ - ra^) = , ml ^ (48) 



For /3 = one has 

m 2 = m i = ? 777^ < , 777q > , 777^ < 777^ < 

or 

m 2 = m 3 = i 777^ < , 777q = , m\{m\ — m\) 7^ 

or 



777, 



= 777^ = , 777^ < , 7774(7774 - 777^) = , 777^ (49) 



We see, that unlike in the previous cases, fine-tuning m\ = alone is never enough to make 
Goldstone sector free of pathologies. So, probably, there is no need to pick out the case 
7772 = i 11 ^ a separate phase. 

3.5 Phase m\ = 

The last remaining option is that 

777^ = , 
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so that vector modes have dispersion relation of the form (26) and do not propagate. This 
case was mentioned in Ref. [15] as a potentially interesting. In this case the field equation 
in the scalar sector still has the form (23) with matrix M being equal to 



2m\v 2 —2%m\v 
1im\v 2{m\ - ml) 

Therefore, the on-shell condition in this case takes the following form, 



d = {4m 2 (ml - ml) - Am\)u 2 = . (51) 

If 

m o( m 3 — m 2) — m 4 = ? (52) 

then matrix M is degenerate, so that one cannot use its inverse as a propagator in the 
perturbation theory and Goldstone sector is not well behaved in the decoupling limit. How- 
ever, for any other values of parameters Goldstone sector is well behaved, and there are no 
propagating vector and scalar modes in the theory (while tensor mode is massive). 

4 Higher- derivative terms, UV (in) sensitivity and un- 
broken gauge symmetries 

In this section we continue the study of the Goldstone sector in the limit when gravity 
is decoupled and discuss how the effects of higher derivative terms affect the conclusions 
of section 3. Let us, however, first discuss what happens when the fine-tuning conditions 
discussed above are slightly violated. Indeed, if these fine-tuning condtions are not ensured 
by symmetries, then they are not stable at the quantum level, so one should check that 
the low-energy effective theory remains well behaved if one allows for small violations of 
fine-tuning relations. The result of such violation is that the non-propagating degrees of 
freedom become dynamical. Then, there are two different situations depending on whether 
the non-propagating degree of freedom had dispersion relation of the form (26) or (27). 

Let us first discuss what happens upon small violation of the dispersion relation (27). 
In the vicinity of the point where this dispersion relation holds, the corresponding degree of 
freedom has dispersion relation of the following form 

uj 2 = ep 2 , (53) 
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where e is a small parameter determined by the distance to the fine-tuned point in the mass 
space. Recall, that the analysis of subsection 3.1 implies that in massive gravity either e < 
leading to classical instability at all momenta, or e > and the corresponding degree of 
freedom is a ghost 6 . However, in both cases these instabilities do no present a problem from 
the viewpoint of the low-energy effective field theory, provided e is small enough. Indeed, 
for negative e the rate T for the development of classical instability is smaller than 

r~cu< |e| 1/2 A 

and can be made arbitrarily small by chosing a small enough |e|, while keeping the cutoff 
value A constant. 

If e is positive, then the corresponding degree of freedom is a ghost, and there is a 
quantum instability related to the possibility to produce it from the vacuum together with 
normal particles of positive energy. However, at small e the maximum (negative) energy of 
the ghost quanta available in the regime of validity of the low-energy effective field theory is 
very small so that the rate of development of quantum instability is strongly suppressed by 
the phase volume. 

Definitely, in both cases there is an instability in the low-energy effective theory, so 
rigorously speaking the latter is not well defined. However, this is not important for any 
practical purposes provided the time-scale for the development of instability is long enough. 
In other words, in all practical situations there is an IR cutoff at long time scales (at least 
at the age of the Universe). 

One could be confused by the fact, that we allow to have ghosts in the low-energy effective 
theory. Recall, (see, e.g., Ref. [18, 19] for a recent discussion) that ghosts imply quantum 
instability of the vacuum due to the possibility of their creation accompanied by normal 
particles of positive energy. In Lorentz-invariant theory the rate of the development of this 
instability diverges due to the infinite phase volume, unless one introduces Lorentz- violating 
cutoff. Then the rate is proportional to the value of the cutoff. In particular, this implies 
that if the only dimensionfull parameter in the theory is the cutoff A itself (as is the case 
in the Goldstone sector, where perturbativity requires that the actual cutoff should be just 
slightly lower than A), then the rate of instability is also of order A. One way to make this 
situation tractable is to introduce a cutoff at the significantly lower scale. 

However, in the Lorentz non-invariant situation it may happen that the dispersion relation 

6 Actually, our analysis does not exclude also a situation, when e is negative and the corresponding degree 
of freedom is a ghost, i.e. both types of instabilities are present simultaneously. 



19 



is such, that the highest (negative) available energy of the ghost quanta is very small (see, 
e.g. dispersion relation (95) which we discuss later). It is even possible, that the energy has 
a maximum at a certain value of momentum and goes to zero at higher momenta. In this 
case processes with substantial production of energy involve a large number of ghosts and 
can be very strongly suppressed even at relatively large values of the coupling constants. It 
may be even possible that there is no need for the explicit Lorentz-violating cutoff to make 
the rate of instability finite — Lorentz-violating dynamics of the system itself may play a 
role of such a cutoff. A somewhat degenerate example of a situation of this type happens in 
quantum cosmology, where scale factor of the Universe is a ghost with zero energy. 

We believe that these issues deserve further studies. In any case, it is clear that if the 
maximum frequency of the ghost quanta is smaller than the inverse age of the Universe, then 
time uncertainty principle makes it nonsensible to speak about production of ghosts. 

Let us now discuss what happens when the dispersion relation of the form (26) is slightly 
violated. In this case one has a mode with dispersion relation 



which is, in our case, again either indicates the presence of classical instability at e < or a 
ghost at e > 0. Here the low values of momenta, 



are dangerous. Indeed, at higher momenta the presence or absence of instability depends on 
the structure of the theory above the cutoff energy, i.e. out of the domain of applicability of 
the low-energy effective field theory, so one can consistently assume that instability is absent. 
To eliminate the dangerously rapid instability at the small values of the momentum one needs 
some IR regulator. This regulator can be provided, for instance, by the finite size of the 
system, or by the Hubble expansion (in some similarity to how Jeans instability is eliminated 
at large length scales). Another option (which, as we discuss in section 5, is realized, e.g., 
for wig = phase) is that mixing with gravity leads to the following modification of the 
dispersion relation (54), 



where m 2 should be positive. 

To summarize, there is an interesting complementarity of two situations (53) and (54). 
Namely, in the first case IR cutoff in time and UV cutoff in space are needed to stabilize the 



2 2 
tU) =p , 



(54) 



P<H 1/2 A 



2 2,2 

euo = p + m , 



(55) 
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theory, while in the second case one needs IR cutoff in space and UV cutoff in time 7 . 

Having discussed what happens upon small departure from the dispersion relations (26) 
and (27) at the one- derivative level, let us turn to the effects of the higher derivative terms. 
Generically, they will also make the would be non-dynamical degrees of freedom dynamical. 
For instance, generic higher derivative terms would lead to the following modification of the 
dispersion relation (27) 



where a is a coefficient of order one. This is precisely the situation in the ghost condensate 
model, where such modification does not lead to any problems (provided a is positive), 
because the resulting degree of freedom has a healthy (positive) kinetic term. However, in 
many of other cases discussed in section 3 this is not the case and the resulting degree of 
freedom is a ghost. Effect of higher derivative terms on the dispersion relation (26) is even 
worse. Indeed, one obtains dispersion relation of the form 



inevitably leading to both quantum and classical instability. 

The above discussion implies that the generic effect of higher derivative terms is to make 
the actual cutoff of the theory much lower than the naive value A. It is worth stressing, 
that this does not imply that one cannot assume that the structure of the theory in the UV 
is such, that coeffcients a, b and their higher order analogues are very small, so that they 
do not lead to any pathologies in the low-energy effective theory below cutoff A, due to the 
same reasons as small departures from fine-tuning relations at the quadratic level. Clearly, 
from the viewpoint of low-energy effective field theory this is a very fine-tuned situation, but 
there is nothing principally wrong in considering such a situation. 

While, a priori, there is nothing wrong in considering a fine-tuned situation, it is worth 
looking for some symmetry principle ensuring the desired fine-tunings and thus making 
massive gravity well behaved without assuming any fine-tuning relations in the Goldstone 
action (besides those, following from the symmetries). An example of massive gravity having 
this attractive property of the UV insensitivity is a ghost condensate. This is achieved by 
requiring that there is a residual invariance under time- dependent spatial reparametrizations 



7 We are thankful to Riccardo Rattazzi for numerous fruitful discussions of these stability issues and for 
the collaboration on the subject. 



2 V 
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(56) 




(57) 
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Besides the purely esthetic attractiveness this property of UV insensitivity may be very 
important from the practical viewpoint as well. Indeed, eventually, we are interested in 
the study of the non-trivial curved backgrounds in massive gravity, and in the absence of a 
symmetry principle there is no guarantee, that, say, the cosmological expansion would not 
drive a sysetm out of the safe region where fine-tuning relations hold. 

To illustrate this, let us consider the following simple choice of function F in the Goldstone 
action (14), 

F = P(X) + (SijYV) 2 (58) 

This is just a ghost condensate action coupled to three free scalar fields 0\ It is straightfor- 
ward to check that fine-tuning condition wig = implies that 

P x + 2XP XX = , (59) 

and that for a large class of functions P one can get a massive gravity in phase wig = 
using this type of actions for Goldstones. On the other hand, equation for field 0° in the 
expanding Universe in this model is the same as in the case of ghost condensate, 

8 t {a*P x d t <f) = , (60) 

implying that at late stages of the expansion a system is driven either to a point 

d t ctP = , (61) 

or to the ghost condensate point 

Px = . (62) 

(in both cases Y goes to zero). However, it is straightforward to check that it is not possible 
to move from the vicinity of the point where condition (59) is satisfied to points (61), (62) 
without crossing regions where instabilities are not suppressed by any small parameter. 

So, let us try to see whether it is possible to protect fine-tuning relations obtained in 
section 3 by residual reparametrization symmetries. In each of the phases found there one of 
the conditions m\ = 0, ml = or m\ = m\ holds. Let us see what symmetries can protect 
each of these relations. 

Let us start from the relation ml = 0. It can be protected by one of the following residual 
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reparametrization symmetries 

x i -> x i + C(t,x) (63) 

x i -> x l + f (i) (64) 

t->* + f°(t,a:) (65) 

t^t + £°(x) (66) 

First symmetry (63) is just a residual symmetry of the ghost condensate, so we will not 
discuss it here. 

At the linear level the second symmetry (64) implies invariance under arbitrary time- 
dependent shifts of the ir 1 fields. In other words, field configuration 

7T° = 0, 7T* = 7r*(t) 

is a solution of linearized field equations, provided there is a residual reparametrization 
symmetry (64). In particular, this implies, that for arbitrary choice of higher derivative 
terms, compatible with symmetry (64), the dispersion relation in the vector sector is (at 
non-zero m|) 

p 2 (l + a^ + 6g + ...)=0, (67) 

where a and b are again some coefficients of order one. This is in agreement with what 
we found at the one derivative level, and implies that there are no propagating light vector 
degrees of freedom in this case. 

Similarly, in the scalar sector symmetry (64) implies, that for generic higher derivative 
terms the on-shell condition in the scalar sector has the following form 

pV» + a^ + . ..) = <>, (68) 

where the lowest order term is determined by our explicit analysis in subsection 3.5. There- 
fore, similarly to the case of ghost condensate, there is one propagating degree of freedom in 
the scalar sector with peculiar dispersion relation 

u 2 oc p 4 (69) 

in this case. In order this degree of freedom had a healthy sign before kinetic term the 
derivative of the eigenvalue of matrix (50) having zero at v 2 = should be positive at this 
point. It is straightforward to check that this is the case, provided 

4 

TYl 
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The point m\ = ml is rather special — here both eigenvalues of matrix (50) have square 
root branchings at v 2 = 0. Perhaps, this case deserves further study. 

Thus, insisting on the residual symmetry (64), one obtains a generalization of the ghost 
condensate theory sharing all attractive features of the latter. On the other hand this model 
has a number of new features. For instance, tensor gravitational waves are massive here. 
Clearly, this theory deserves a separate detailed study. 

Symmetry (65) implies that all masses but ml and m\ vanish. As we discussed in sub- 
section 3.5 in this case matrix M is degenerate, so that Goldstone sector is not well behaved 
in the decoupling limit. 

Finally, symmetry (66) does not allow to put any restrictions on higher-derivative correc- 
tions to the dispersion relation in the vector sector, so it does not allow to avoid fine-tuning. 
Similarly, in the scalar sector it implies that dispersion relation has a form 

w V + a— + ...) = , 

that does not protect from the instabilities. 

Let us now discuss symmetries protecting fine-tuning relation m\ = m\. Besides sym- 
metry (63), leading to ghost condensate, this fine-tuning relation can be protected by a less 
restrictive residual symmetry 

+ (71) 
Actually, this symmetry implies a stronger condition 

mj = ml = m \ = . (72) 

In the vector sector it protects from any instabilities in the low-energy effective theory, 
leading to the dispersion relation 

c 2 (l + a^ + ...), (73) 

provided m\ ^ 0. In the scalar sector symmetry (71) implies the following form of the 
dispersion relation (at the lowest order we are using the result of subsection 3.3) 

cV + a^ + ...) = 0. (74) 

Consequently, we again obtain one propagating degree of freedom in the scalar sector with 
dispersion relation (69). To understand when this mode has a healthy kinetic term let us 
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a) 



b) 



Figure 1: Splitting of zeroes of the eigenvalue of the matrix M due to higher-derivative 
terms in the presence of residual symmetry (71) in two different cases: a) m\ , Wq > 0, b) 

m\ , ml < 

look at the behavior of the eigenvalues of matrix (24) when masses m%, m\ and m\ are zero. 
At this point one of the eigenvalues has a double zero at uo 2 = while another has no zeros at 
all. Higher order terms in the dispersion relation (74) lead to the splitting of the double zero 
into two simple zeroes as shown in Fig. 1 (for negative value of the coefficient a in Eq. (74)). 
Now, for positive m\ the propagating mode has a healthy kinetic term in the situation shown 
in the left panel of Fig. 1, while for negative m\ in the situation shown in the right panel of 
Fig. 1. Straightforward calculation shows that propagating mode is not a ghost provided 



This UV insensitive theory resembles ghost condensate even more, because tensor mode is 
massless in both cases. Also, in both cases graviton mass has a form of a gauge fixing term 
and thus higher order terms (non-linear, or higher derivative) are needed to see modification 
of gravity. 

Finally, let us discuss whether it is possible to make use of a residual gauge symmetry 
to avoid fine-tunings in the phase with ml=0. Besides the symmetry (65) leading to the 
degeneracy in the Goldstone sector one may try the following symmetry 



m\m\ > . 



(75) 



(76) 
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This symmetry implies that 

m o — m 4 — • 

Unfortunately, in the scalar sector this symmetry implies the dispersion relation of the form 

j9V + a^ + ...) = (77) 

and thus is not enough to obtain a UV insensitive theory. 

To summarize, we found two new examples of the UV insensitive massive gravities. This 
is the whole phase mf = and a part of the boundary of m\ = m\ phase, where condition 
(72) holds. In both cases the residual symmetry group is a subgroup of the residual symmetry 
of the ghost condensate. Both theories deserve further detailed study. 

Clearly, our analysis does not exhaust all possible subgroups of the diffeomorphism group. 
It is worth studying whether other interesting possibilities exist. 



5 Mixing with gravity and absence of the vDVZ dis- 
continuity 

Let us discuss now, how mixing with gravity may affect the above conclusions and show that 
the vDVZ discontinuity is absent if graviton masses are in the ranges found in section 3. 

Let us start from the discussion of the vDVZ discontinuity. The complete linearized field 
equations in massive gravity can be presented schematically as follows 

(M ° +M ~i h :;) = { T °)' (78) 

where T^ u is the energy-momentum tensor of matter, we assumed that there are no sources 
for the Goldstone fields ir^; matrix _M is diagonal 

/ Ml,V% \ 
*°=( ADS J ' (?9) 

and matrix M. m accounts for graviton masses and kinetic mixings between Goldstones and 
gravitons, 

/ A 4 A 4 V m , 
" ~ A 4 £> m I' <80) 
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Here we assumed that gauge symmetry is fixed in the graviton sector (for instance, using 
harmonic gauge), so that T>\ is the ordinary two-derivative graviton kinetic operator fol- 
lowing from the Einstein-Hilbert action with gauge fixing terms, T>q is the two derivative 
kinetic operator in the Goldstone sector, O is a dimensionless matrix, determining the ten- 
sor structure of the graviton mass matrix, and T>, m is a one-derivative operator describing 
mixing between Goldstones and graviton. For graviton mass parameters found in section 3 
both operators T>\ and T>q are invertible at general values of energy 00 and momentum p. 
Therefore, one can present matrix M. in the following form 

^(^(ISS A ^)K <8i) 

This presentation makes it manifest, that, far from the on-shell values of energy and mo- 
mentum found in section 3, the metric induced by the source T^ v differs from that in the 
pure Einstein theory by small corrections, vanishing when the scale of graviton masses 

A 2 

m 



M Pl 

goes to zero. In particular, this is true in the Newtonian regime, provided energy uj is much 
larger than the graviton mass scale 

m < uj < p < A . 

Hence, on time and length scales shorter than the inverse graviton mass one reproduces the 
usual Newtonian limit of gravity and vDVZ discontinuity is absent. 

Let us check now whether mixing with gravity can introduce dangerous instabilities for 
the "safe" phases found in section 3. The most straightforward way to do this is, following 
Ref. [15], to study the spectrum of massive gravity in the unitary gauge in each of the cases. 
Namely, one writes 

h 00 = i> (82) 
hoi = Ui + diV (83) 
hij = Xij + (diSj + djSi) + didjcr + d^r (84) 

where ^ is a transverse-traceless tensor, Ui and Vi are transverse vectors and other fields are 
three-dimensional scalars. Let us start from the vector sector. Here one has the following 
action [15] 

L v = Mp t (sid^djSi - UidjUi + 2w i <9 2 <9 Sj + m\uiUi + m 2 2 s i d 2 j s l ) . (85) 
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Dispersion relation following from this action is 

2 2 2 2 2/2, 2\ n 

m-tf) 00 — m 2 p yp + m x ) = (J . 

Consequently, the only effect of mixing with gravity in the vector sector is the emergence of 
the mass gap m\. 

Let us now consider the scalar sector. Here, at general values of masses, one has the 
following quadratic action in the unitary gauge 

M 2 

L = (4 (if> - 2d v + dla) d 2 r + 6rd 2 r - 2rd 2 r + mfy 2 + 2m 2 1 (d i v) 2 - 

™\ ( 2 + 2rSf <t + 3r 2 ) + m\ (d 2 a + 3r) 2 - 2mfy (d 2 a + 3r) ) . (86) 

Let us start our analysis without assuming any fine-tuning relations between masses. Fields 
if) and v enter action (5) without time derivatives [15] (the term propartional to dovdfr in 
Eq. (5) may be written as vd dfr). Field equations resulting from variation with respect to 
these fields can be written as 

v = 4,<V (87) 
m\ 

a=- 2 T-^T + ^ 2 U. (88) 
mi Of mi of 

Plugging expressions (87) and (88) back into action we arrive at the following Lagrangian for 
the two remaining fields r and ip (we will discuss later what happens if m\ = or m\ = 0), 

* - ^ ((£ - A) * * - ^ - 6 * + (M - 2 ) f T - fe ^ + 

i^Tdli, - i'4 (ml - ml) rdU' - Ul + ^(m? - mi)) i? + im^A ■ < 89 > 
m% ml x \ ml J mi J 

Let us first consider the phase m\ = 0, when all terms in the second line of Eq. (89) vanish. 
But before actually setting them to zero, the following comment is in order. Straightforward 
inspection of action (89) shows, that the dispersion relation of the mode which is non- 
dynamical at wig = takes the following form at small (much smaller than all other masses) 
no n- vanishing m 2 ,, 

mlu 2 = -|^(p 2 + ^ , (90) 



where parameter /Iq is defined as follows 

4 



2 2 2 

mi m{ 



28 



We see, that mixing with gravity automatically provides an IR regulator of the type (55) 
needed to stabilize the Goldstone sector upon small departure from the surface m\ = 0. 

Stability of the action (89) in the bulk of the m\ = phase was studied in Ref. [15]. As 
we mentioned in subsection 3.2 the only difference of the results of this analysis with ours is 
that extra stability condition 

Ami > m l ( 91 ) 

was imposed in that work. If this condition is violated, the term proportinal to rdfr enters 
action (89) with negative sign. Depending on other parameters, this may lead to the Jeans- 
like instability in the finite region of small momenta and frequencies, u,p ~ m. Presence 
of such instability does not invalidate the use of the low-energy effective theory with cutoff 
A for time scales shorter than ~ m" 1 and may even lead to interesting phenomenological 
consequences. So we do not think that one should a priori disregard mass values violating 
condition (91). 

Let us see now what are the effects of mixing with gravity at the boundaries of the phase 
mo = 0. If m\ = 777.4, then action (89) describes a propagating degree of freedom with the 
following dispersion relation (cf. Ref. [15]) 

u 2 = (p 2 + zfij)^ + ml, (92) 

where 

^7779— 777q 3 ,777n „ „ . . 

2 2 4 3 = -, 4-f-l = 3z. 93 

ml /if 777| 

We see, that in this case mixing with gravity invalidates our conclusion that there are no 
propagating degrees of freedom in the scalar sector. The new degree of freedom has a healthy 
kinetic term and does not lead to the vDVZ discontinuity provided m\ 7^ m\. In order to 
avoid rapid classical instabilities, one needs, however, 

or, equivalently, 

7772 > 7773 . (94) 

This condition was impossible to obtain from inspecting the Goldstone Lagrangian in the 
decoupling limit at m\ = m\. 

If 7772 = m i the dispersion relation obtained in Ref. [15] reduces to 

- = ^^TW ■ (95) 
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Again we see that mixing with gravity results in the emergence of a new propagating degree 
of freedom. This degree of freedom does not lead to any rapid instabilities, provided 

nl > , 

or, equivalently, 

m\ > mj > . (96) 

As before, this condition is impossible to deduce from the Goldstone Lagrangian alone. 
Finally, if in addition to ml = one has ml = 0, then Eq. (88) implies that r = 0. Then 
action in the scalar sector takes the following form 

L = ^f{m\-m\){d}af (97) 

and does not describe any propagating degrees of freedom. Consequently, mixing with gravity 
does not bring nothing new in this case. 

Let us now discuss effects of mixing with gravity for phase m\ = 0, protected by the 
residual gauge symmetry (63). Here Eq. (87) implies, that either uj = 0, or r = 0. In the 
latter case, one uses Eq. (88) to eliminate tp and arrives at the following action for a, 




which does not describe any propagating fields. So in this case mixing with gravity does not 
introduce any extra degrees of freedom. 

We believe, that the above discussion is representative enough to demonstrate what are 
the possible effects of mixing with gravity. So we leave the detailed study of the consequences 
of this mixing in the m\ = m\ phase and for m\ = with negative m\ beyond the scope 
of the current paper. Instead, we just consider the point corresponding to the second UV 
insensitive theory found in section 4. 

Before doing that, let us mention one more possible effect of mixing, that was not realized 
in the phase m\ = 0. Namely, it follows from Eq. (89) that in the phase with m\ = ml (and 
non-zero ml) a generic consequence of mixing with gravity is that sixth polarization mode of 
the graviton becomes dynamical. Discussion in subsection 3.3 suggests that this mode may 
be a ghost or lead to classical instabilities. Presumably, there are regions of mass parameters 
where this is not dangerous. Our discussion in section 4 implies, that what one needs, is 
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that all quanta of this mode are very soft (cf. dispersion relation (95)) and, consequently, 
the corresponding instability rate is very slow 8 . 

For point m\ = m| = m| = 0, protected by the symmetry (71), one uses Eq. (88) to 
express ip through r and Eq. (87) to express v through r. As a result one arrives at the 
following action 

L = ^k (±. ( d f T y _ A(9 ^r) 2 + Ad^rd^a + 6(d r) 2 + 2(d t r) 2 ) . (99) 

Z, y /77q 777. ^ J 

Variation of this action with respect to a implies that r = if both uj and p are non-zero. 
Then variation with respect to r gives a — 0. So mixing with gravity does not introduce 
new propagating modes in this case as well. 

6 Discussion 

In this final section let us first briefly discuss what kind of cosmology can one expect in 
models of massive gravity discussed above. A hope that long scale modification of gravity 
may explain the observed acceleration of the cosmological expansion rate, or, at least, relate 
it to effects testable at shorter scales (such as the anomalous precession of the Moon [20, 21]) 
was the main driving force for the interest to massive gravity and its brane world analogues 
during the last few years. The common lore, supported by the brane world models (e.g., 
by the DGP model [10]) is that modification of the gravitational potential at the distance 
scale r c is related to the modification of the Friedmann equation at the value of the Hubble 
expansion rate H of order r^ 1 . 

An example of the ghost condensate suggests that in four- dimensional models of massive 
gravity situation may be completely different. Indeed, at the linear level this model is char- 
acterized by distance scale r c and time scale t c , such that gravitational potential of static 
sources is strongly modified at distance and time scales larger than r c and t c correspondingly. 
On the other hand, the homogeneous Friedmann equation in this model is completely equiv- 
alent to the standard one (modulo, possibly, the presence of additonal dust-like component, 
with relative density depending on the initial conditions for the scalar field) for arbutrary 
value of the Hubble rate. 

To illustrate that this is not a peculiar property of the ghost condensate, let us consider 

8 Note, that this may imply additional limitations on the possible values of masses, that are not seen in 
the decoupling limit. 
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Lorentz- violating gravity with function F (see Eq. (14)) of the following form 

F = A 4 F{XY 3 ) . 



(100) 



In the unitary gauge this means, that we are adding to the Einstein action the following 
non-covariant term 

J rf 4 xA 4 v^F(^°(Tr^) 3 ) . 

We will comment later what was the reason to assume this particular form of the action. 
At the moment, let us just note that independence of function F on Z can be ensured, e.g., 
by time inversion symmetry. Also, for the solutions we discuss here V 1 = 0, so possible 
dependence of F on V 1 would not affect our discussion. Finally, there is a global symmetry 

protecting dependence of function F on the combination XY 3 at the quantum level. 

Flat homogeneous cosmologies in massive gravity are described by the following ansatz 
for the metric and Goldstones, 

ds 2 = dt 2 - a 2 (t)dx i dx i (101) 
0° = 4>°(t) (102) 
f = AV (103) 

In particular, Y = 1 /a 2 and goes to zero as the Universe expands, so one may suspect that 
effect of all masses but wiq is negligible at late times. We will see, however, that this is 
not necessarily the case. Ansatz (101) automatically satisfies (f) 1 field equations, while 0° 
equation reads as follows 

d t (a 3 F'(XY 3 )Y 3 d t (j) ) = dt (F'(Xa- 6 )X 1/2 a~ 3 ) = . (104) 
A solution to this equation is 

X = x a 6 

or 

0o (t) = x l J 2 J dta 3 (t) . 

Combination XY 3 entering the Goldstone action remains constant for this solution. On the 
first hand this implies that graviton masses are constant during the cosmological evolution. 
In particular, if some fine-tuning relation were true in the begining of expansion it keeps to 
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be true during the whole history of the Universe. For instance, it is straightforward to check, 
that action (100) allows cosmological solutions in the = phase of massive gravity. 

On the other hand, time-independence of XY 3 implies that contribution of the Goldstone 
sector to the energy- momentum tensor is also constant. In other words, this contribution is 
either zero or cosmological constant. So this is an example of the model where all graviton 
modes are massive, but flat homogeneous solutions are the same as in the Einstein gravity. 
Should one be disappointed by this fact? We believe that the answer is not. This just 
implies that the value of the graviton mass may be quite large, affecting growth of primordial 
perturbations, structure formation, or even gravitational dynamics at the (super)galactic 
scales. 

Even more intriguingly, this opens up a possibility that the effects of backreaction of 
small scale inhomgeneities on the expansion rate of the Universe can be quite substantial in 
massive gravity. Actually, the backreaction effects in conventional Einstein theory are not 
completely understood, but it appeares that the likely conclusion is that they are very small 
in the realistic Universe (see, e.g. Refs. [22, 23, 24, 25] for a recent discussion). However, 
this may be completely different in massive gravity. Indeed, in this case there is at least one 
extra length scale in the problem — inverse mass of the graviton. Now imagine the Universe 
consisting of the heavy point-like galaxies at large distances from each other. Then, when 
the distances between neighboring galaxies are larger than the inverse graviton mass, the 
pairwise interactions between galaxies are practically negligible (at least, in the rriQ = 
phase, where the gravitational potential is of the Yukawa type [15]). Is it possible to use in 
this case the homogeneous approximation in cosmology? The answer to this question is not 
clear to us. 

An indirect hint that there may be a subtlety here is coming from the attempt to find 
closed or open homogeneous cosmological solutions in the model (100). In this case instead 
of a distance between galaxies there is another extra length scale in the problem — curvature 
radius 9 . It is straightforward to see, however, that the corresponding ansatz does not satisfy 
all the equations (except the case of the de Sitter space, when there are extra symmetries). 
This implies that closed or open Universe is necessarily inhomogeneous on the scale of its 
curvature radius in massive gravity with action (100). 

To conclude this brief discussion of the cosmology in massive gravity, it is worth noting 
that time-independence of the graviton mass during the cosmological expansion is by no 
means a general property of massive gravity. Actually, the particular form of the action (100) 
9 We thank Dominik Schwarz for suggesting us to look for open or closed solutions. 
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was derived from the requirement that graviton mass is constant. The study of cosmological 
consequences of the time-dependent graviton mass may provide new surprises. 

Finally, let us stress that there is a large number of problems not addressed in the current 
paper. First of all, there are still a lot of questions about stability of massive gravity. For 
instance, our analysis of mixing with gravity in section 5 covers only the phases wig = 
(where it only slightly extends that performed in Ref. [15]), m\ = and a point with 
enhanced symmetry in the phase with m\ = m\. More detailed analysis of the phases with 
m 2 = m i an d m 2 = with negative m\ is needed to decide which regions of parameters 
identified in section 3 actually correspond to tractable low-energy effective theories with 
cutoff A. Also, it is necessary to include higher derivative terms in this analysis. A particulary 
interesting question is whether two new UV insensitive theories found here may avoid Jeans- 
like instability present in the ghost condensate case. Furthermore, the analysis of stability 
of non-linear solutions (on the first hand, cosmological ones and massive analogues of the 
Schwartzchild metric) is needed to conclusively decide on the phenomenological acceptability 
of models discussed here. 

Also, in the current paper we concentrated on the issues of consistency and stability and 
leave aside the analysis of phenomenology of modification of gravity in models studied here. 
Definitely, one of the most pressing questions is what are the actual limits on the graviton 
mass parameters. It seems likely, that to obtain an accurate answer to this question a detailed 
study of growth of primordial perturbations and structure formation in massive gravity is 
needed. 

Another circle of questions is what are the limits on the direct couplings between Gold- 
stone sector and Standard Model fields and what are the implications of massive gravity for 
inflation (see, Refs. [26, 27] where these issues were studied in the case of ghost condensate). 
We hope to address these and other related questions in future. 
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